/CR-1998-206848 


<3 & % 3 3 6 ' 




SeaWiFS Science Algorithm Flow Chart 


Michael Darzi 



National Aeronautics and 
Space Administration 

Goddard Space Flight Center 
Greenbelt, Maryland 2077 1 


January 1998 






The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key 
part in helping NASA maintain this important 
role. 

The NASA STI Program Office is operated by 
Langley Research Center, the lead center for 
NASA’s scientific and technical information. The 
NASA STI Program Office provides access to 
the NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA’s institutional 
mechanism for disseminating the results of its 
research and development activities. These 
results are published by NASA in the NASA STI 
Report Series, which includes the following 
report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data or 
theoretical analysis. Includes compilations of 
significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA’s counterpart of 
peer-reviewed formal professional papers but 
has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, techni- 
cal, or historical information from NASA 
programs, projects, and mission, often con- 
cerned with subjects having substantial public 
interest. 

• TECHNICAL TRANSLATION. 
English-language translations of foreign scien- 
tific and technical material pertinent to NASA’s 
mission. 

Specialized services that complement the STI 
Program Office’s diverse offerings include creat- 
ing custom thesauri, building customized data- 
bases, organizing and publishing research results . 

. . even providing videos. 

For more information about the NASA STI Pro- 
gram Office, see the following: 

• Access the NASA STI Program Home Page at 
http://www.sti.nasa.gov/STI-homepage.html 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access Help 
Desk at (301) 621-0134. 

• Telephone the NASA Access Help Desk at 
(301) 621-0390. 

• Write to: 

NASA Access Help Desk 
NASA Center for AeroSpace Information 
800 Elkridge Landing Road 
Linthicum Heights, MD 21090-2934 






NASA/CR-1998-206848 


SeaWiFS Science Algorithm Flow Chart 


Michael Darzi, SAIC General Sciences Corporation, Laurel, Maryland 


CONTRACT NAS5-32377 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 2077 1 


January 1998 



Available from: 


NASA Center for AeroSpace Information 
800 Elkridge Landing Road 
Linthicum Heights, MD 21090-2934 
Price Code: A17 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 
Price Code: A10 



introduction 


This flow chart describes the baseline science algorithms for the SeaWiFS Data Processing 
System (SDPS). As such, it includes only processing steps used in the generation of the 
operational products that are archived by NASA's Goddard Space Flight Center (GSFC) 
Distributed Active Archive Center (DAAC). It is meant to provide the reader with a basic 
understanding of the scientific algorithm steps applied to SeaWiFS data. It does not include 
non-science steps, such as format conversions, and places the greatest emphasis on the 
geophysical calculations of the Level-2 processing. Finally, the flow chart reflects the logic 
sequences and the conditional tests of the software so that it may be used to evaluate the 
fidelity of the implementation of the scientific algorithm. In many cases however, the chart 
may deviate from the details of the software implementation so as to simplify the 
presentation. 

The flow chart was compiled as part of the SeaWiFS Project Calibration/Validation 
Element's effort to evaluate the scientific algorithms used for processing SeaWiFS data. 
The following people also contributed to the compilation of the flow chart: Wayne D. 
Robinson, Eueng-nan Yeh, Robert E. Eplee, Norman A. Kuring, and Bryan A. Franz. 

A World Wide Web version of this chart is available as part of SeaWiFS's main web site. 
The web version is updated whenever algorithm changes are implemented. The chart was 
designed for web use with hyperlinks connecting its pages. 

Flow Chart Legend 

A rectangular box represents a process, which, when shaded, indicates that an expanded 
chart of that process is depicted in another figure. A diamond shaped box indicates a 
binary decision and always has two exit paths representing "yes" and "no" results. A circle 
or oval represent a Boolean combination of paths. The DAAC is depicted by an octagon 
representing the repository for all operational products. Finally, a shaded triangular box 
with a lowercase letter indicates that a comment exists relevant to the processing or logic 
depicted in that area of a figure. The letter identifies the comment within the set of 
comments given for the figures of the flow chart. 

A small dashed line represents the functional limit of each figure's topic. Straight lines 
indicate paths of science information, such as primary data. Dashed straight lines indicate 
optional or occasional science information paths. Curved lines indicate control paths, 
without science information transfers. 

Italicized text describe data being transferred along the straight line paths or, in larger font, 
the starting source data. Normal, small font text is used to provide explanations of paths 
when needed. 
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Comments 

The following comments accompany the figures of the flow chart. Each set of comments 

are labelled with the title of the figure to which they apply and the page number of that 

figure in the chart. 

SeaWiFS Data Processing (2): 

a. Only global area coverage (GAC) data are processed operationally beyond Level-1 . 

b. Each GAC scene is normally one swath of data taken on the sunny side 

(descending node) of an orbit. : 

c. The space binning step converts each swath into binned products. 

d. The time binning step combines (i) all space-binned scenes of a day into a day- 

binned file; (ii) each 8-day set of consecutive day-binned files (starting each year 
from Julian day 1) into an 8-day binned product; (iii) all the day-binned products of 
each calendar month into a month-binned product; and (iv) all month-binned 
products of each calendar year into a year-binned product. Thus, day- and month- 
binned products are time binned again to generate the longer period time-binned 
products. 

Ancillary Data Conversion (3): 

a. For each SeaWiFS scene, the ancillary data closest in time to before and after the 
scene, and during the scene, if any, are selected for processing with that scene. 
Each ancillary data parameter is interpolated to the location and time of each pixel 
in the scene. (See L2 Support Data Calculations, p.11.) If in the unusual 
situation where such near-real time ancillary data are not available, climatologies of 
the ancillary parameters are substituted. 

b. Interactive examination occurs as part of normal quality control procedures or when 
an automatic statistical check program, through which all ancillary data are run, 
indicates possibly problematic data. 

c. Bad data imply unrealistic or missing data. Correction of such data involves an 
analyst replacing the problematic areas with data from the climatologies or other 
averages, or using gridding or smoothing techniques. 

L0-L1A Conversion (4): 

a. Recorded data consist of GAC and local area coverage (LAC) data. Each GAC 
scene is normally one swath of data taken on the sunny side (descending node) of 
an orbit. Each LAC scene, including calibration scenes, normally consists of one 
continuous recording of high-resolution data. For HRPT data, each Level-0 
collection forms one scene. 
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LI Browse Generation (5): 

a. The calibration applied is the same as that shown under Sensor Corrections (p.7) 
without the stray-light correction. 

b. The navigation of pixels is performed using data stored in the Level-1 A product 
during L0-L1A Conversion (p.4). 

c. The 8-bit image stored in a Level-1 browse file can be converted to a 24-bit (8 bits 
red, 8 bits green, 8 bits blue) image by application of a color look-up table. The 
8-bit values stored in the browse file are no more than indices into that look-up table 
and should not be construed to bear any direct relationship to radiances measured 
by the sensor. Once the look-up table has been applied, the red, green, and blue 
components of a pixel will be close to--but probably not the same as-the scaled 
Rayleigh reflectances (for the 670, 555, and 412 nm bands, respectively) computed 
by the browse file generator before it quantized the 16,777,216 possible colors down 
to 256 or fewer. 

L2 Processing (6): 

a. For LAC resolution (not archived products for Levels-2 and -3), 194 anchor points 
are defined. 

b. A flag may be designated as a "mask" by program input. 

c. The final assignments of the 12 geophysical parameter values that are output in the 
Level-2 products are shown on the following pages: 

- page 15: La_865 and eps_78; 

- page 18: nLw_412, nLw_443, nLw_490, nLw_510, nLw_555, and La_670; 

- page 20: K_490; 

- page 23: tau_865; 

- page 24: CZCS_pigment and chlor_a. 

d. The alternate values for the 12 geophysical parameters are: 


Geophysical Parameter 

Alternate Value 

nLw_412 

Level- 1 A radiance counts of band 

nLw_443 

Level- 1 A radiance counts of band 

nLw_490 

Level-1 A radiance counts of band 

nLw_510 

Level- 1 A radiance counts of band 

nLw 555 

Level- 1 A radiance counts of band 

La 670 

Level- 1 A radiance counts of band 

La_865 

Level- 1 A radiance counts of band 

CZCS_pigment 

0 

chlor a 

0 

K_490 

0 

eps_78 

0 

tau_865 

0 
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A summary of the conditions assigned to the various Level-2 flags, and the pages 
of the flow chart on which the assignments occur, is given in the following table: 



Flag 

Page 

QpngjtjQn 



1 

13 

invalid tilt state 



1 

13 

band-4 Rayleigh less than or equal to 0 

I 

| 


1 

15 

aerosol determination error 

i 

j 

1 

15 

epsilon out of range 

i 

j 

1 

18 

La(1 ..6) less than or equal to 0 

■ 

1 

j 

2 

13 

land 

j 

1 

3 

13 

problematic ancillary data 

j 

i 

4 

13 

ZGLINT greater than threshold (glint) 

: 

I 

5 

7 

L1A(1..8) greater than knee value 

= 

: 

6 

13 

satellite zenith greater than threshold 


5 

7 

13 

shallow water 


1 

8 

13 

missing bands or bad navigation 

J 

] 

8 

13 

L1A(1..8) less than or equal to 0 


j 

8 

18 

Lw(1..6) less than or equal to 0 

1 

f 

8 

18 

nLw(1 ..6) less than or equal to 0 

1 

5 

9 

8 

stray light 

- 

| 

10 

13 

ice or cloud 

| 


10 

24 

not shallow water & flag 8 is set 


■■ 

11 

19 

coccolithophores J ~- 

S 


12 

22 

T6 greater than T5 (turbid water) 

= 


13 

13 

solar zenith greater than threshold 


_ 

14 

23 

tau_865 greater than threshold 

= 

— 

15 

18 

nLw(5) less than threshold 


’ 

16 

24 

chlorophyll algorithm error 



Sensor Corrections (7): 

a. The calibration table is occasionally updated to incorporate new information about 
the sensor's performance. The table always contains the calibration information for 
all data collected since the start of the mission. 

b. A time-dependent gain and an offset may be specified as program input to override 
those read from the calibration table. This is normally done for testing purposes 
only. 
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f^'^onawct'stray’tah. correction is requested by program input (normal option), 
the GAC stray-light correction routine will work with a rotating buffer of 3 scansu 

hrat k (^test) C s^ t 'ln n me\uff^w r hwe i r ) can ^, 1 ^ ' Hluen^ by me' center 
products), a similar process is used with a 5-scan line buffe . 

“ •"^StSS? JSmfemtoSe anchors (defined on p.6. L2 Processing) on both 
sides of, and on the same scan as, the current pixel. 

1 2 Calculations (A) (13): . 

a. A flag may be designated as a "mask" by program input. 

i 2 Calculations (B) (14): . . 

a. A flag may be designated as a "mask" by program mpu . 

Atmospheric Correction (15): 

a. A flag may be designated as a "mask" by program input. 

Radiance Calculations (18): 

a. A flag may be designated as a "mask" by program input. 

Coccollthophore Test (19): 

a. A flag may be designated as a "mask" by program input. 

a' 9men Note 1 that the PGMT value calculated here Is only used for “I® L' U -Xn a*!* CZCS 
tho next page. The pigment value calculated on p.2 , P y . ^ 

Pigment , T the one used as a geophysical parameter output in the Uevel-2 

products. 

Turbid Water Test (22): . 

a. A flag may be designated as a "mask" by program input. 

Aerosol Optical Thickness (23): 

a. A flag may be designated as a "mask' by program input. 
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Chlorophyll a & CZCS Pigment (24): 

a. A flag may be designated as a "mask" by program input. 

Space Binning (26): 

a. The 13 geophysical parameters that are space binned are the 12 parameters of 
Level-2 products (see comment (c) under L2 Processing for a list) plus 
chlor_a_K_490, calculated during space binning. 

b. The mean of any geophysical parameter in a space-binned product is equal to 
MEAN bJ = SUMX bJ / W b for bin b and parameter j. The standard deviation is equal 
to {[(SUMXX b /W b )-MEAN bJ 2 ]W b 2 /(W b 2 -S b )} 05 for bin b and parameter j. Note that 
space-binned products are not archived. 

Time Binning (27): ■ ~ 

a. The time binning step combines (i) all space-binned scenes of a day into a day- 

binned file; (ii) each 8-day set of consecutive day-binned files (starting each year 
from Julian day 1) into an 8-day binned product; (iii) all the day-binned products of 
each calendar month into a month-binned product; and (iv) all month-binned 
products of each calendar year into a year-binned product. Thus, day- and month- 
binned products are time binned again to generate the longer period time-binned 
products. - • - ~ 

b. The 13 geophysical parameters that are time binned are those of the space-binned 
products (see comment (a) under Space Binning). 

c. The mean of any geophysical parameter in a time-binned product is equal to 
MEAN bJ = SUMX bJ / W b for bin b and parameter j. The standard deviation is equal 
to {[(SUMXX b /W b j-MEAN b /]W b 2 /(W b 2 -S b )} 05 for bin b and parameter j. 

Standard Mapped Image Generation (28): 

a. A different standard mapped image product is generated for each of five 
geophysical parameters, chlor_a, CZCS_pigment, nLw_555, tau_865, and K_490, 
from each time-binned product. 

Scaling for Byte Values (30): 

a. The 8-bit image stored in a Level-1 browse file can be converted to a 24-bit (8 bits 
red, 8 bits green, 8 bits blue) image by application of a color look-up table. The 
8-bit values stored in the browse file are no more than indices into that look-up table 
and should not be construed to bear any direct relationship to radiances measured 
by the sensor. Once the look-up table has been applied, the red, green, and blue 
components of a pixel will be close to-but probably not the same as-the scaled 
Rayleigh reflectances (for the 670, 555, and 412 nm bands, respectively) computed 
by the browse file generator before it quantized the 16,777,216 possible colors down 
to 256 or fewer. 
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Ancillary Data Conversion 




Source Ancillary Data 




^ EPTOMS 03 

1 

NCEP 

A 

TOVS 03 ^ 

1 

Winds, Pressure, & 
Predpitable Water (PW) 

1 
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L0-L1 A Conversion 


LO Data 
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L2 Processing 


03, Winds, 
Pressure, A PW 


— ry pat 
uonversig 


L1A Data 


(GAC) 




If any flag 1,8,10,11,12,14, 
1 5,1 6 is set & is a "mask" 


Use alternate values 
(radiance counts or 0) 
for all L 2 fields 



I 


Space Binninc 

, 



L2 Data 
(GAC) I 
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Stray-Light Correction 
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Stray-Light Detection 



Test for a BT 


if [DELTA > LTF r (L1 B(8,P+1) - Ltyp)] and [LI B(8,P+1) > 0.9*knee_value], 
then a left edge (LE) has been detected at P+1 of S; 


if [-DELTA > LTF*(L1 B(8,P) - Ltyp)] and [LI B(8,P) > 0.9*knee_value], 
then a right edge (RE) has been detected at P of S; 
where LTF is a fractional value specified by program input (baseline value is 
0.25), Ltyp is a typical sea-surface radiance for band 8 (1 .09 mW/(cm2 um sr)) f 
and where P and P+1 are not already both marked as BT. 
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Bright-Target Processing 



The left edge (LE) of a BT has 
been detected at P+1 of S. 


Left edges 


The right edge (RE) of a BT has 
been detected at P of S. 


Mark up to 3 pixels (that are not marked as BT or AT) 
to the left of LE by their distance from LE. 



Keep Incrementing P and testing for RE until RE is found or there are 
no more pixels left In S: 

1) if LI B(8,P) < 1 .25*Ltyp or 

2) If [-DELTA > LTF*(L1 B(8,P) - Ltyp)] and [L1B(8,P) > 0.9*knee_value], 
then a right edge (RE) has been detected at P of S. 



Mark up to 3 pixels (that are not marked as BT or AT) 
to the right of RE by their distance from RE. 



If no LE had been found since the start of S or since the previous RE, 
keep decrementing P and testing for LE until LE is found or there are 
no more pixels left in S: 

1) f LI B(8,P) < 1 .25*Ltyp or 

2) f [DELTA > LTF*(L1 B(8,P+1) - Ltyp)] and 
[L1B(8,P+1) > 0.9*knee_value], 

:hen a left edge (LE) has been detected at P+1 of S. 


If associated LE Is not found, assume that first pixel in scan is the LE. 


Assume that the last pixel In scan is the RE. 
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L2 Support Data Calculations 



03, Winds, 
Pressure, & PW 


LI A Data & 
LIB Data 



Calculate day-of-year correction, FOCORR, to solar constant 
and the solar constant , FO, for all 8 bands 


Interpolate pixel latitudes & longitudes (from 
L0-L1 A Conversion) to anchor points 


Interpolate pixel satellite zenith angles (from 
L0-L1 A Conversion) to anchor points 


Calculate the solar zenith angles for the anchor points 


Calculate the cosine of the solar zenith angles for the anchor points 


For ancillary data (03, surface U & V winds, surface pressure, & 
precipitable water), do 2D & time interpolation for each pixel in scan 
then interpolate those pixel values to values for the anchor points 


ANQUAL = 0 for all anchor points. 

If an anchor point’s ancillary data were derived with data flagged as problematic 
during Ancillary Data Conversion, ANQUAL <> 0 for that anchor point 


At eacn ancnor point, calculate windspee 
> WS = fUwind**2 + Vwind**21**0.5 


Use pressure and viewing geometry to calculate 
Rayleigh reflectances for the anchor points 


Use WS and viewing geometry to calculate glint 
coefficients as fraction of F0(8) for the anchor points 


Calculate normalized water-leaving whitecap reflectance at each anchor point: 
> WHITE = [6.49*1 0**(-7)] * WS**3.52 


Use pressure and satellite geometry to calculate the Rayleigh component 
of the satellite/surface absorption for the anchor points 


Use pressure and sun geometry to calculate the Rayleigh component 
of the sun/surface absorption for the anchor points 


Use 03 and viewing geometry to calculate satellite/surface/sun 
ozone correction for the anchor points 


Use 03, pressure, and viewing geometry to calculate satellite/surface/sun 
ozone and Rayleigh correction for the anchor points 


current scan indicates a new sensor 
tilt state or if the tilt is changing, calculate 
Fresnel reflectance at the anchor points 
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L2 Support Data Interpolations 



JBSEiSmal 




ZLAT = pixel latitude interpolated from adjacent anchors | A 

ZLON = pixel longitude interpolated from adjacent anchors] 


ZTHETA = pixel satellite zenith angle interpolated from ad 


ZTHETO = pixel solar zenith angle interpolated from adjacent ancnorsi 

ZMUO = pixel cosine of solar zenith angle interpolated from adjacent anchors 


ZANQUAL = 0 

If ANQUAL of either anchor adjacent to the pixel <> 0, ZANQUAL <> 0 


ZRAYLY = pixel Rayleigh reflectance interpolated from ad 


ZGLINT = pixel glint coefficient interpolated from adjacent anchors 


ZWHITE = WHITE of adjacent anchors interpolated to pixel 


ZB1 * Rayleigh component of the satellite/surface absorption for pixel, 
interpolated from adjacent anchors 


ZB2 = Rayleigh component of the sun/surface absorption for pixel, 
interpolated from adjacent anchors 


ZBS03 * satellite/surface/sun ozone correction 

for pixel, interpolated from adjacent anchors 


ZBST = satellite/surface/sun ozone and Rayleigh correction 
for pixel, interpolated from adjacent anchors 


ZFRESNL = Dixel Fresnel reflectance interpolated from adjacent anchors 
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L2 Calculations (A) 


Apply whltecap correction to total radiances (LIB data), Lt, of bands 1 to 8 
> Lt(B) = Lt(BL- ZWH IT E*ZBST (B)*F0( B)*ZM UQ/PI where B = 1 ..8 


Apply 03 correction to Lt for bands 1 to 8: 
> Lt(B) = Lt(B)*ZBS03 where B = 1 ..8 


^^Are ZLON 
ZLAT within the land 
''s. mask? 


If ZLON & ZLAT are within shallow water mask, set flag 7 


ZMUCT 
v> 0?> 


Band-8 albedo = 100*Lt(8)*ZB1(8)*ZB2(8VF0(8 


If albedo > Input threshold (Implying ice or clouds), set flag 10 


If ZTHETO (solar zenith angle) > Input threshold, set flag 13 


If ZTHETA (satellite zenith angle) > input threshold, set flaa 6 


If scan-line quality flag #2 (from L0-L1 A Conversion) Indicates 
a missing band, or if bad-navigation flag is on. set flaa 8 


If It Is not a valid tilt state (from L0-L1A Conversion), set flag 1 


If ZANQUAL (problematic ancillary data) <> 0, set flag 3 


If any band-4 Rayleigh reflectance at the adjacent anchor 
points In the current or previous scan Is <= to 0. set flaa 1 


If LI A count for any band Is <= to 0, set flag 8 


If ZGLINT > Input threshold, set flag 4 (sun glint 




If any flag 1 > 2.3.4,5.6.7.8.9, 
10,13 la set & la a "mask* 


Use alternate values 
(radiance counts or 0) 
for all L2 fields 
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L2 Calculations (B) 
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Atmospheric Correction 



If It Is the first pixel of a new scan or If the last successful aerosol determination 
was for a pixel more than 12 (for GAC) or 48 (for LAC) pixels away from the 


current one, force a new search for bounding aerosol models by making the 
previous model indices "unavailable" (see Aerosol Determination) 


MU0 = cosine of ZTHETO 


Convert Lt to reflectances and subtract Rayleigh reflectances for bands 1 to 8: 
> REFL(B) = { PI*Lt(B)/[F0(B)*MU01) - ZRAYLY(B) where B = 1 ..8 


Apply 02 correction to REFL(7 


Aerosol Determination 


Get the stored diffuse transmittance coefficients, a & b, for the bounding models. 


Compute the diffuse transmittance for a bounding set 
of view angles using TAU from Aerosol Determination: 
> t(B) = a * exp[-b °TAU(B)1 where B=1..8 


Use interpolation ratio from Aerosol Determination 
to calculate t at the pixel 


For bands 1 to 6, calculate water-leaving radiances times transmittance: 
> tLw( B) = (REFL(B) - RHO_A(B)1*FO(B)*MUO/PI where B = 1..6 


, calculate aerosol radiances: 

> La(B) = REFL(B)*F0(B)*MU0/P1 where B = 7,8 

(La(8) is used as geophys. parameter La_865 In the L2 products) 

If an error has occurred during fZZZZZZZZUZZI^” 

Aerosol Determination, set tlagll ~ 

Using the interpolation ratio and the derived I if f| ag f j S set 

epsilon of the bounding models from Aerosol if flag 1 1s set \ & is a ’mask' 
Determination, calculate epsilon (geophys. & is a 'mask' ] 
parameter eps_78 In L2 products) and, if N. / 

outside Input threshold range, set flag 1 \ /A 



Use alternate values 
(radiance counts or 0) 
for all L2 fields 
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Aerosol Determination 


Atmospheric Correctlo 


Compute the aerosol single & multiple scattering relationship coefficients for all 


nal an error & skip rest of aerosol determination 


I Kg:U] ftlcffl BlUO EJ lajigll 1 k 



/''Are the lndices v \ 
tor the bounding model! 

used for the previous 
“\nixel available?/' 


Get the single scattering aerosol terms for 
the 2 bounding models at bands 7 & 8 using 
the scattering relationship coefficients and 

preclpltable water (PW) 

if undefined operation occurs, signal error I 
and skip rest of aerosol determination 

Calculate the epsilons for the 2 bounding 
models using the band 7 to 8 ratio of these 
scattering terms and average them 


Aerosol Bounding Mo dels 


/ is the average’'^ 
within the EPS values 
of the 2 bounding,^ 
models? 


Calculate the aerosol optical thickness, TAU, for the 8 bands and the 2 bounding 
models, using the view geometry and the band-8 single scattering reflectances 


Calculate the Interpolation ratio of the average 

epsilon from those of the 2 bounding models | 

Use the epsilons of the 2 bounding models, the single scattering term for band 8, 
the PW, the scattering relationship coeffs., and the Interpolation ratio to derive 
the multiple scattering aerosol reflectance terms, RHO_A, for bands 1 to 6 


Use the single scattering aerosol terms for the 2 bounding 
models at bands 7 & 8 and the interpolation ratio to compute 
the actual single scattering terms for bands 7 & 8 
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Aerosol Bounding Models Search 


M = 1 



et single scattering reflectances for model M at bands 7 & 8 i 
RHO_A(7) & RHO_A(8), the scattering relationship coeffs, and 


If undefined operation occurs, signal error 
and skip rest of aerosol determination 


Calculate the epsilon for the band 7 to 8 ratio of 


the single scattering reflectances for model M 


Add model M epsilon to sum of epsilons for models<M 


Average the derived epsilons (divide sum 


Calculate the differences between the average 
epsilon and each of the 12 band 7/8 model epsilons 


N = 10 


Rank the differences according to their absolute values 


Sum the N smallest differences 


to the N smallest differences and weighted 
In an Inverse proportion to those differences 


Average sum of the derived epsilons 


Calculate the differences between the average 
epsilon and each of the 12 band 7/8 model epsilons 


N = N - 2 


8 using 
ind PW 


:m > 12?; 


Select as the lower bound model the one whose difference with 
the 4 derived epsilon average Is the smallest positive value 

Select as the upper bound model the one whose difference with 
the 4 derived epsilon average Is the smallest negative value 


Retain Indices of bounding models for the next pixel 
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Radiance Calculations 


AtmosphericCorrection 


Calculate band-6 aerosol radiance: 

> La(6) = Lt(6) - [F0(6)*ZMU0*ZRAYLY(6)/Pi] - tLw(6) 
(used as geophys. parameter La_670 in L2 products; 
see L2 Calculations (A) for Lt calculation) 


Calculate water-leaving radiances, Lw, for bands 1 to 6 (remove transmittance): 
> Lw(B) - tLw(B) / t(B) where B = 1..6 


!MU0 > 0’ 


inLw(B) = 0 where B = 1 ..6 


Calculate normalized water-leaving radiances, nLw, for bands 1 to 6: 

> nLw(B) = Lw( B)*ZB2( B)/(ZM UO* FOCORR) where B = 1 ..6 

(geophys. parameters nLw_412, nLw_443, nLw_490, nLw_510, nLw_555 

in L2 products) . 


If Lw(B) <= 0 for B= 1 ..6, set flag 8| 

If nLw(B) <= 0 for B= 1 .. 6 , set flag 8 j 
If La(B) <= 0 for B = 6..8, set flag 1 f— 
< input threshold, set flag 1 5 


\WSM 


If any flag is set 
that is a 'mask* 


If flag 15 is set 
& is a "mask* 



Use alternate values 
(radiance counts or 0) 
for all L2 fields 
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Coccolithophore Test 



SeaWiFS Science Algorithm Flow Chart 


Science information flow Control path Functional limit of topic 


January 5, 1998 Page 1 9 






iililWHWIIIIUIIi Hill 


Diffuse Attenuation 
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Pigment 


Diffuse Attenuation 


^ Are nLw(2), ' 

nLw(3), and nLw(5) 
\ all >0? 


PGMT = 0| 


PGMT = ex 


PGMT >yes_ 


PGMT = 


[nLw(3)/nLw( 


0.592 - 3.48*[nLw(3)/nLw(5)] 


JO < PGMT < 2? 


■nLw(2)/n Lw(5) -^no 
^j> 0.163? ^ 


PGMT = 0| 


PGMT = PGMT*1 .280*sart 


PGMT^yes. 

■v<0?/ 




0.163 



Turbid Water Test 
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Turbid Water Test 
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Aerosol Optical Thickness 
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Chlorophyll a & CZCS Pigment 


Aerosol Optical Thickness 


RRS(5)<0 or 
*RRS(3)<0?. 


Convert nLw for bands 1 to 6 to reflectances: 

> RRS(B) = nLw(B)*FOCORR/FO(B) for B = 1..6, 

M chlor_a =~0 

PGMT = 0 




Calculate L2 products’ geophys. parameter chlor_a = 

-0.40 + 10**(0.341 - 3.001 *RATIO + 2.81 1 *RATIO**2 - 2.041 RATIO 3) 


2hlor_a<=0.01?, 

'■'Tno 


chlor_a = 0.01 


PGMT = 1.34 1 

k 0.01**0.98 


PGMT = 1 .34 * chlor_a**0.98 

(used as geophys. parameter CZCS_plgment in L2 products) 


ilor_a > input threshold?" 
^current threshold^-- 


7 is not set (not shallow water) and 
is set fnLw<=0l set flag 10 (cloud/ice) 


Set flag 1 6 


If flag 10 is set 
& is a ’mask* 


If flag 16 is set 
& is a ’mask* 


L2 Processing 


Use alternate values 
(radiance counts or 0) 
for all L 2 fields 
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Space Binning 



Check if scene crosses data day boundary. 

■ « _ • _ l I J.. !«■ 


If so, use separate bin grids (data days) for 
binning eastern and western portions of scene. 
I (Two binned scene products will be output.) 


I Initialize sums to 0| 


Get next scan 


Get next pixel 


V es ^More\ n ° 
\scans?X I 


More 

Dixels?^ 


X Are any L2 or 
,X engineering quality \ 
flags set that are designated 
v for exclusion or is a bad > 
^.navigation or bad tilt X 
\ indicated? 


Y no (bin 
this pixel) 


Calculate bin index for pixel’s lat/lon. 
For that bin index: 

> increment pixel count (N) 

> set scene count (S) to 1 

> chlor_a_K_490 = chlor_a/K_490 

> For each of 13 geophys. parameter 

> add pixel value to sum (SUMX) 

> add pixel value A 2 to > 

sum of values A 2 (SUMXX) ‘ 


yes (skip 
this pixel) 


For each bin with N>0: 

> set time trend (TT) to 1 

> calculate the weight (W) 
as the square root of N 

> For all geophys. params: 

> SUMX = SUMX/W 

> SUMXX = SUMXX/W 


L3 Binned 

Data 

(scenes) 
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Time Binning 
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Scaling for Byte Values 



m Browse 

B = byte value to form true-color Image (L 1) 


B = [log(D)+2] / (2/255) 
ml B = 0 <==> D = 0.01 

/& B = 255 <==> D = 1 .0 

or to write Into product (L2 & L3) 

D = Rayleigh-corrected reflectances (L 1) 
or geophysical value (L2 & L3) 


H2 Browse Generation! 

B = [log(D)+2] / 0.01 5 for chlor_a, and if B>250, B=250 
B = 0 <==> D = 0.01 mg 

B = 250 <==> D = 56.234133 mg m^3 

{Standard Mapped image (SMI) Generation! 

B = [log(D)+2] / 0.01 5 for chlor_a & CZCS_pigment, and if B>254, B=254 

B = 0 <==> D = 0.01 mg m A -3 

B = 254 <==> D = 64.565423 mg m^S 

B = D / 0.063 for nLw_255, and if B>254, B=254 

B = 0 <==> D = 0.01 mW cm A -2 um'M si^-l 

B = 254 <==> D = 16.002 mW 017^2 um A -1 sr A -1 

B = D / 0.005 for tau_865, and if B>254, B=254 

B = 0 <==> D = 0.01 

B = 254 <==> D = 1 .270 

B = D / 0.025 for K_490, and if B>254, B=254 
B = 0 <==> D = 0.01 m*-1 

B = 254 <==> D = 6.350 nY'-l 

[ L3 Browse Generation ! 

B = SMI chlor_a byte value and if 250<B<255, B=250 
B = 0 <==> D = 0.01 mg 0^3 

B = 250 <==> D = 56.2341 33 mg rrA-3 


Summary o f Reserved Bvta Values 


Byte 

LI brs 

L2brs 

SMI 

1 3 hrs 

- 

0-250 

data 

data 

data 

data 


251 

data 

not used 

data 

caption (1)* 


252 

data 

flag 4 (3) 

data 

geocoordinate gnd (3) 


253 

data 

flag 2(1) 

data 

land (4)* 


254 

data 

flags 1.5, or 10 (2) 

data 

political boundaries (2)* 

■ 

255 

data 

not used 

no data 

no data (5) 

■ 


n Numbers Indicate "priority of alignment: lowest vaiue => highest "priority I 
* Indicates feature that Is exercised on/y Ifre^uestedb^mgra^n^^^ 
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